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The general solution of the Einstein equation for higher dimensional (HD) spherically symmetric 
collapse of inhomogeneous dust in presence of a cosmological term, i.e., exact interior solutions of 
the Einstein field equations is presented for the HD Tolman-Bondi metrics imbedded in a de Sitter 
background. The solution is then matched to exterior HD Scwarschild-de Sitter. A brief discussion 
on the causal structure singularities and horizons is provided. It turns out that the collapse proceed 
in the same way as in the Minkowski background, i.e., the strong curvature naked singularities form 
and that the higher dimensions seem to favor black holes rather than naked singularities. 
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I. INTRODUCTION 

That the results comingfrom the analysis of high red- 
shift Type la supernovae |3, 0, B indicate the Universe is 
accelerating. This suggests the possibility that a nonzero 
cosmological constant (A) may dominate the total energy 
of our Universe. The cosmological implications of the ex- 
istence of a A term are enormous, concerning not only the 
evolution of the Universe, but also structure formation 
and age problems. If A term must be restored to the Ein- 
stein equations, surprises may turn up in other physical 
applications of Einstein's equations as well. For example, 
Markovic and Shapiro 0| generalized the Oppenheimer- 
Snyder model (which describes the gravitational collapse 
of a spherical homogeneous dust ball initially at rest in 
exterior vacuum to a Schwarzschild black hole) taking 
into account the presence of a positive A. They showed 
that A may affect the onset of collapse and decelerate 
the implosion initially. The results of the Markovic and 
Shapiro were qualitatively generalized to the inhomoge- 
neous dust and degenerate cases by Lake for both 
A > and A < 0. It was shown explicitly by Cissoko et 
al. 1^ that the cosmological term slows down the collapse 
of matter, limiting the size of black holes. 

The Tolman-Bondi metric 7| , which is asymptotically 
flat, has been extensively used to study the formation of 
naked singularities in spherical collapse. It is seen that 
the Tolman-Bondi metric admits both naked and covered 
singularities depending upon the choice of initial data 
and there is a smooth transition from one phase to the 
other However^ccording to 

the cosmic censorship conjecture (CCC) [131, the singu- 
larities that appear in gravitational collapse are always 
surrounded by an event horizon. Moreover, according 
to the strong version of the CCC, such singularities are 
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not even locally naked, i.e., no non-spacelike curve can 
emerge from such singularities (see 18], for reviews on 
the CCC). The study in the inhomogeneous dust collapse 
with a positive A, from the viewpoint of C CC. was ex- 
amined in [a, 123, l21| ■ Deshingkar et al. [l^l showed 
that the presence of a A can cover a part of the sin- 
gularity spectrum which is visible in the corresponding 
dust collapse models for the same initial data. Whereas 
Wagh and Maharaj [22| showed that in spherically sym- 
metric radiation collapse (Vaidya collapse), the effect of 
adding a positive A does not radically alter the descrip- 
tion. Lemos arrived at the same conclusion for a 
negative A. The result in both these cases are the same 
as in the case of collapsing radiation in the Minkowskian 
background. Therefore, at least, in the case of spherical 
radiation collapse, the asymptotic flatness is not essential 
for the development of a naked singularity. 

While gravitational collapse has been originally stud- 
ied in four dimensions (4D), there have been several at- 
tempts, mainly motivated by strin g th eory, to study it in 
HD space-time [lllMElliiESlliilHSEll- Since, 
current experimental results involving tests of the inverse 
square law do not rule out extra dimensions even as large 
as a tenth of a millimeter. It is now important to consider 
the evolution of the extra dimensions since the observed 
strength of the gravitational force is directly dependent 
on the size of the extra dimensions. As a consequence, 
there is a renewed interest towards understanding of the 
general relativity in more than four dimensions, as grow- 
ing volume of recent literature indicates. In particular, 
several solutions to the Einstein equations of localized 
sources in higher dimensions have been obtained in the 
recent years j32^. 

In this paper, we shall study spherical inhomogeneous 
dust collapse with a positive A in HD theory of gravity, 
and present solutions in closed form. This is HD anal- 
ogous of 4D Tolman-Bondi-de Sitter solutions and for 
definiteness we shall call it HD Tolman-Bondi-de Sitter 
solutions. Then, we show that HD Tolman-Bondi-de Sit- 
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ter admits strong curvature naked singularity. However, 
the presence of a positive A does not radically alter the 
established picture of Inhomogeneous dust collapse. 

In the next section, we give exact HD spherically sym- 
metric solution of Einstein field equation for a collaps- 
ing inhomogeneous dust with a cosmological constant A. 
This is followed by junction conditions between a static 
and a non-static HD spherically symmetric space-time in 
section HI. The nature of singularities of such a space- 
time, and the consequence of cosmological constant A is a 
subject of section V. This is preceded by detailed analysis 
on apparent horizon in section IV. 

We have used units which fix the speed of light and 
the gravitational constant via SttG = — 1. 



II. HIGHER DIMENSIONAL TOLMAN-BONDI 
DE SITTER SPACE-TIMES 
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where the overdot and the prime denote the partial 
derivative with respect to t and r, respectively. Inte- 
gration of Eq. (O gives 
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which can be substituted into Eq. © to yield 
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The standard 4D Tolman-Bondi solution represents 
an interior of a collapsing inhomogeneous dust sphere. 
The solution we seek is - collapse of a spherical dust with 
a positive A in HD space-time. We choose a spherically 
symmetric comoving metric in HD |23l |25| . which has 
form 



ds^ = df - e^^^'^^^dr^ - R{t,rfdn^ 
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is the metric on an n-sphere and n = D ~ 2 (where D 
is the total number of dimensions), together with the 
stress-energy tensor for dust: 



Tab = Cit,r)Sl5l 



(3) 



where Ua = is the (n -I- 2)-dimensional velocity. The 
coordinate r is the co-moving radial coordinate, t is the 
proper time of freely falling shells, and i? is a function of 
t and r with R > and A is also a function of t and r. 
With the metric the Einstein equations are 
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The functions AA{r) and /(r) are arbitrary and referred 
to as the mass and energy functions, respectively. Since 
in the present discussion we are concerned with gravita- 
tional collapse, we require that R{t, r) < 0. The energy 
density C,{t, r) is calculated as 



C(tr) = 



nM' 
2i?"i?' ' 



(10) 



For physical reasons, one assumes that the energy den- 
sity C(ij f) is everywhere non-negative. The special case 
/(r) =0 corresponds to the marginally bound case which 
is of interest to us in this paper. Substituting Eqs. (jHJ 
and 10 into Eq. yields 



M' = -Ci?"i?'. 



Integrating Eq. leads to 
M{r) = ■ 



(11) 
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where constant of integration is taken as zero since we 
want a finite distribution of matter at the origin r — Q. 
The function Ai{r) must be positive, because M.{r) < 
implies the existence of negative mass. This can be seen 
from the mass function M(i,r), which is given by 
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The quantity Ai{r) can be interpreted as energy due to 
the energy density ({t, r) given by Eq. ((T^ . and since it 
is measured in a comoving frame, M is only r dependent. 
Cissoko et al. have derived marginally bound (/ = 0) 
dust solution in the presence of A > 0. Here we derive the 
analogous HD solutions. Equation jnj, for vanishing A, in 
4D as well as in HD, has three types of solutions, namely, 
hyperbolic, parabolic and elliptic solutions depending on 
whether /(r) > 0, f{r) = or f{r) < 0, respectively. 
The condition /(r) = and A = is the marginally 
bound condition, limiting the situations where the shell is 
bounded from those it is unbounded. In the presence of a 
cosmological constant, the situation is more complex, and 
/(r) = leads to an unbounded shell. The assumption 
allows for analytical solutions in closed form 
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where a = a{t, r) has the form 
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where io(^) is an arbitrary function of integration which 
represents the proper time for the complete collapse of a 
shell with coordinate r. It follows from the above that 
there is a space-time singularity at i? = and at i?' = 0. 

It is easy to see that as A ^ the above solution 
reduces to the HD Tolman-Bondi solutions 1251: 
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The standard 4D Tolmon-Bondi solution can be now re- 
covered by setting n — 2. The three arbitrary functions 
A4{r), f{r) and to{r) completely specify the behavior of 
shells with radius r. It is possible to make an arbitrary 
relabeling of spherical dust shells by r g{r), without 
loss of generality, we fix the labeling by requiring that, 
on the hypersurface t = 0, r coincides with the radius 

i?(0,r) = r. (20) 

This corresponds to the following choice of to{r): 
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The central singularity occurs at r = 0, the correspond- 
ing time being t = to{0) = 0. We denote by p{r) the 
initial density: 

n M' 2 /■ 

p{r) = C(0,r) = M{r) = - p{r)r^dr. (22) 



III. JUNCTION CONDITIONS 

In order to study the gravitational collapse of a finite 
spherical body we have to match the solution along the 
time like surface at some R — i?s to a suitable HD ex- 
terior. We consider a spherical surface with its motion 
described by a time- like (n-l-l)-surface S, which divides 
space-times into interior and exterior manifolds V/ and 
Ve- According to the generalized Birkoff theorem 
the vacuum space-time outside is HD Schwarzschild-de 
Sitter space-time: 

ds^ = F{Y)dT^ - J^dY^ - y^dn^, (23) 

where is a function of Y given by 

2M 2Ay2 
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and M is a constant. In accordance with Darmois junc- 
tion condition, we have to demand when approaching E 
in Vi and Ve 



(25) 



where the subscript S means that the quantities are to 
be evaluated on S and let Kf^ is extrinsic curvature to 
S, defined by 



dx 



(26) 



and where Tp^ are Christoffel symbols, the unit nor- 
mal vectors to E, x" are the coordinates of the inte- 
rior and exterior space-time and arc the coordinates 
that defines S. The intrinsic metric on the hypersurface 
r = rs is given by 



ds^ = dt^ - i?2(rs,i)dr2^ 



(27) 



with coordinates — {t, 6i, 02, 03, . . , On)- In this co- 
ordinate the surface E, being the boundary of the matter 
distribution, will have the equation 



0, 



(28) 



where rs is a constant. The first fundamental form of E 
can be written as gijd^'^d^^ . Then the exterior metric, on 
E, becomes: 
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where we assume that the coefBcient of dT^ > so that T 
is time hke coordinate. From the first junction condition 
we obtain 



1 (dY^\ 
F{Y^) [dT J 
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The non- vanishing components of extrinsic curvature 
of S can be calculated and the result is 
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With the help of Eqs. (jSOI) - l|SlIl and © , the total 
energy entrapped within the surface S can be given by 

M 



n — 1 
M = ^—M{r). 



(35) 



Thus, the junction conditions demand that the HD 
Schwarzschild mass M is given by Eq. . 



IV. HORIZONS 

The apparent horizon is formed when the boundary of 
trapped n spheres are formed. In spherical dust collapse, 
the event horizon coincides with the apparent horizon at 
the boundary of the spherical mass distribution. The 
apparent horizon is the solution of 



g'^'R^aR.b - -R' + f{r) + 1 
Upon using Eqs. |(SJ) and ©, we have 
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For A = we have the Schwarzschild horizon i?"^^ = 
A4, and for = we have the de Sitter horizon 
R — ±y/n{n + l)/2 A. The approximate solutions of 
Eq. H37|l to first order are 
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where Rhh is the radius of the black hole event horizon. 
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where Rch is the radius of the cosmological event horizon 
|35| . There exist a critical solution of Eq. (|37|l . where two 
roots coincides and there is only one horizon. The time 
for the formation of apparent horizon, from Eq. I|15|) is 
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Equation (|41|l . gives the time for the formation of event 
horizon in HD Tolman-Bondi-de Sitter space-time. A 
necessary condition for the singularity to be globally 
naked is Iah > to{0). 



V. CAUSAL STRUCTURE OF SINGULARITIES 

It has been shown ^| that Shell-crossing singularities 
are characterized by i?' = and R > 0. On the other 
hand the singularity at i? = is, where all matter shells 
collapses to a zero physical radius and hence known as 
shell focussing singularity. 

We shall consider the case t > to. In the context of 
the Tolman-Bondi models the shell crossings are defined 
to be surfaces on which R' = {R > 0) and where the 
density C diverges. A regular extrcmum in R along con- 
stant time slices may occur without causing a shell cross- 
ing, provided C{t,r) does not diverge. By Eq. (^UJ, this 
implies A4' = where ever R' = and also that the 
surface R' = remain at fixed R. Now Eq. Hl()|) implies 
to = 0. Thus the condition for a regular maximum in 
R{t, r) is that M' = 0, t'^ = hold at the same R. It 
has been shown j^^l that shell crossing singularities are 
gravitationally weak and hence such singularities cannot 
be considered seriously. 

Next, we turn our attention to shell-focusing singular- 
ities. Christodoulou pointed out in the 4D case that 
the non-central singularities are not naked. Hence, we 
shall confine our discussion to the central shell focusing 
singularity. The energy density diverge at t = tQ{r) in- 
dicating the presence of curvature singularity 37]. It is 
known that, depending upon the inhomogeneity factor, 
the 4D Tolman-Bondi solutions admits a central shell 
focusing naked singularity in the sense that outgoing 
geodesies emanate from the singularity. Here we wish 
to investigate the similar situation in our HD Tolman- 
Bondi-de Sitter space-time. We consider a class of mod- 
els such that 



7W(r) =7r"-\ 
to{r) = Br. 



(42a) 
(42b) 



This class of models for 4D space-time is discussed in 
[ll[li|. The parameter B gives the inhomogeneity of the 
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collapse. For _B = all shells collapse at the same time. and the equation of general density becomes 

For higher B the outer shells collapse much later than the 

central shell. We are interested in the causal structure 

of the space-time when the central shell collapses to the 

center (i? = 0). From Eqs. ifTUI) and H42al) . the energy 

density at the singularity 
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As t ^o(^) i-e. in approach to singularity, we have 
sinh a w a and coth a w 1/a then Eq. (|44|l reduces to: 
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The nature (a naked singularity or a black hole) of the 
singularity can be characterized by the existence of radial 
null geodesies emerging from the singularity. The singu- 
larity is at least locally naked if there exist such geodesies, 
and if no such geodesies exist, it is a black hole. The 
critical direction is the Cauchy horizon. This is the first 
outgoing null geodesic emanating from r = t = Q. The 
Cauchy horizon of the space-time has y = t/r = const 
[lOl [l5| . The equation for outgoing null geodesies is 



dt 
dr 



= R'. 



Hence along the Cauchy horizon, we have 
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and using Eqs. H47|l and H16() . with our choice of the 
scale, we obtain the following algebraic equation: 
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To facilitate comparison with Ghosh and Beesham |25j, 
we introduce a a new relation between B and 7 as, 



B 
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(49) 



with the help of above relation and after some rearrange- 
ment, Eq. H48() takes the exactly the same form as in 
Ghosh and Beesham 1231 



n + 1 



\/lV 



V72/-l = 0. (50) 



This algebraic equation governs the behavior of the tan- 
gent vector near the singular point. The central shell 
focusing singularity is at least locally naked, if Eq. H5U|) 
admits one or more positive real roots. Hence in the ab- 
sence of positive real roots, the collapse will always lead 
to a black hole. Since the Eq. (|50l) is the same as that 
obtained by Ghosh and Beesham 25] when the metric is 
asymptotically flat i.e., A = 0. Consequently the values 
of roots for the geodetic tangent and the condition for 
these values to be real and positive are the same as those 
obtained for the asymptotically flat situation in '25^1 . As 
a result, the Tolman-Bondi-dc Sitter space-time has same 
singularity behavior as the Tolman-Bondi space-time in 
both 4D and HD. Thus the results of collapsing inho- 
mogeneous dust in de Sitter background are similar to 
that of collapsing inhomogeneous dust in Minkowskian 
background, as it should have been expected, since when 
t to{r) the cosmological term 2AR^/n{n + 1) is negli- 
gible. It is known that the Tolman-Bondi metric (A = 0), 
in the 4D case is extensively used for studying the for- 
mation of naked singularities in spherical gravitational 
collapse. It has been found that Tolman-Bondi metric 
admit both naked singularities and black holes form de- 
pending u pon the choice of initial data. Indeed, both 
analvtical |llj|-[T^ and numerical results 8] in dust in- 
dicate the critical behavior governing the formation of 
black holes or naked singularities. One can now safely 
assert that end state of 4D Tolman-Bondi collapse is now 
completely known in dependence of choice of initial data. 
A similar situation also occurs in HD Tolman-Bondi col- 
lapse [23 125II for a class of model under discussion. Since 
the presence of A > does not change the final fate of In- 
homogeneous dust collapse. Hence, we can conclude both 
naked singularities and black holes can form in the HD 
spherical inhomogeneous dust collapse with a positive A. 

Further, it is shown that, in the class models discussed 
here, the formation of a black hole is facilitated with in- 
troduction of the extra dimensions or, in other words, the 
naked singularity spectrum gets continuously covered in 
higher dimensional space-times. This is valid in the in 
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null fluid collapse |^ and scalar field collapse |28|| in 
higher dimensional space-times. Thus it appears that the 
singularity will be completely covered for very large di- 
mensions of the space-time. Recently, It has been shown 
that cosmic censorship can be restored in some classes of 
inhomogeneous dust gravitational collapse models when 
space-time dimension to be A'^ > 6, i.e., the naked singu- 
larities of this model can be removed by going to higher 
dimensions |30| . However, this is valid with only smooth 
initial profiles. 

VI. CONCLUDING REMARKS 

In this paper. We have shown that the 4D spheri- 
cally symmetric solution describing inhomogeneous dust 
collapse with a cosmological term go over to (n + 2)- 
dimensional spherically symmetric solution and essen- 
tially retaining its physical behavior and when n — 2, 
one recovers the 4D Tolman-Bondi-de Sitter solutions 
ll^ I2H . Thus we have obtained Tolman-Bondi-dc 
Sitter metric in arbitrary dimensions and the junction 
condition for static and non-static space-times are de- 
duced. We have also utilized this solution to study the 
end state of collapsing star and showed that there exists 
a regular initial data which leads to a naked singularity. 

Our purpose was to investigate the collapse of inhomo- 
geneous dust shells in an expanding de Sitter background, 
to find out if the naked singularity occurs in this situa- 
tion and to compare any difference with the similar col- 
lapse in the asymptotically flat case. We have obtained 



a condition for the occurrence of a naked singularity in 
the collapse of dust shells in an expanding background 
which is the same as that obtained when the background 
is asymptotically flat. This very fact establishes that the 
space-time is asymptotically flat or not docs not make 
any difference to the occurrence of a naked singularity. 
This is evident at least the in the class of models defined 
by Eq. (g^l- 

Penrose has conjectured that it seems unlikely that 
a A-term will really make much difference to the singu- 
larity structure in a collapse. The relevance of A is really 
only at the cosmological scale. Our results are consistent 
with this. However, The introduction of a cosmologi- 
cal constant changes scenario in many ways. There are 
now several apparent horizons instead of one. However, 
only two apparent horizon are physical, namely the black 
hole horizon and the cosmological horizon. Other results 
derived in Q do carry over to HD space-time essentially 
with same physical behavior hence not presented to avoid 
duplication. 

Finally, the result obtained would also be relevant in 
the context of superstring theory which is often said to 
be next " theory of everything" , and for an interpretation 
of how critical behaviour depends on the dimensionality 
of the space-time. 
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